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ABSTRACT 

Stellar substructure within the Milky Way’s halo, including dwarf galaxies, star clus¬ 
ters, and streams, holds a wealth of information about the formation and evolution 
of our Galaxy. The detection and characterization of these structures require homoge¬ 
neous deep imaging data over a large area of the sky. The Dark Energy Survey (DES) 
is a 5000 sq. degree survey in the southern hemisphere, which is rapidly reducing the 
existing north-south asymmetry in the census of MW satellites and other stellar sub¬ 
structure. We use the first-year DES data down to previously unprobed photometric 
depths to search for stellar systems in the Galactic halo, therefore complementing the 
previous analysis of the same data carried out by our group earlier this year. Our search 
is based on a matched filter algorithm that produces stellar density maps consistent 
with stellar population models of various ages, metallicities, and distances over the 
survey area. The most conspicuous density peaks in these maps have been identified 
automatically and ranked according to their significance and recurrence for different 
input models. We report the discovery of one additional stellar system besides those 
previously found by several authors using the same first-year DES data. The object is 
compact, and consistent with being dominated by an old and metal-poor population. 
DES J0034-4902 is found at high significance and appears in the DES images as a 
compact concentration of faint blue point sources at ~ 87kpc. Its half-light radius of 
rh = 9.88 ± 4.31 pc and total luminosity of My ~ —3.05Jq; 42 are consistent with it 
being a low mass halo cluster. It is also found to have a very elongated shape. In ad¬ 
dition, our deeper probe of DES 1st year data confirms the recently reported satellite 
galaxy candidate Horologium II as a significant stellar overdensity. We also infer its 
structural properties and compare them to those reported in the literature. 

Key words: Stellar populations; Galaxy evolution; stellar statistics 


2 E. Luque et al. 


1 INTRODUCTION 

The census of Milky Way (MW) satellites has grown rapidly 
over the past fifteen years. Several of these newly found ob¬ 
jects are star systems with very low luminosities (—3.0 < 
Mv < 0) and small half-light radii (< 10pc), being more 


consistent with star clusters l 

Koposov et al.|2007 Belokurov 

et al.| 2010; jFadely et al. 2011[|Munoz et al. 2012; Balbinot 

et al. 2013 

1. These clusters are thought to be suffering stellar 


mass loss via dynamical processes such as tidal disruption 
or evaporation ( Koposov et al.|2007 Kim & Jerjen 2015a). 
The number of dwarf galaxies around the MW has also in¬ 
creased significantly, from the 12 classical dwarfs known un¬ 
til the late 1990s, up to a total of 27 which were known by 
early this year (McConnachie |2012 1, thanks in large part to 
the Sloan Digital Sky Survey (SDSS). Several of these lat¬ 
ter were shown to be very low luminosity systems with high 
M/L, thus representing some of the most dark matter rich 
objects. 

At the larger luminosities typical of globular clusters 
(GCs), different cluster subpopulations classified by their 
position, kinematics and horizontal branch morphology have 
been known for several decades (Zinn 1985 1993). The so- 
called young halo clusters may have originated in dwarf 
galaxies accreted by the MW ( Zinn|[l993 Lee et aJ7|| 1994 ). 
Both types of objects seem to share a vast planar structure 
around the Galaxy, which also encloses several stellar and 


gaseous streams of clusters and dwarf galaxies (Pawlowski 
|et al.|[2012[ |2015| ) . The accretion origin of part of the MW 
system of GCs is also supported by the fact that several of 
them are found to have positions and kinematics that re¬ 
late them to the Sagitarius dwarf galaxy (|Law fe Majewski| 


2010). On the other hand, at the much lower luminosities of 


the recent satellite discoveries, the very distinction between 
star clusters and dwarf satellites may become less clear, as 
attested by their respective loci in size and luminosity space. 
It is therefore important to pursue a complete census of 
faint stellar systems inhabiting the Galactic halo, and to 
characterize them in terms of structure, stellar populations 
and dark matter content. Extrapolations of the SDSS results 
over the entire sky and over the currently known luminosity 
function of MW dwarfs indicate that this census is still very 
incomplete ( Tollerud et al.|2008 Hargis et al.|2014 i. 

A very recent boost to the number of known MW satel¬ 
lites has been brought by the Dark Energy Survey (DES; 
The Dark Energy Survey Collaboration! [2005). Using the 


first infernal release of DES coadd data (Y1A1), Bechtol 


et al. (20151 reported on the discovery of eight new MW 


satellites over a solid angle of 1,800 square degrees in the 
southern equatorial hemisphere. Six of these systems have 
sizes and optical luminosities clearly consistent with the low- 
luminosity dwarfs previously detected in SDSS. The case for 
the other two objects is less clear. In a parallel effort, |Ko-| 


posov et al. (20151 reported nine new MW satellites using 


the same DES imaging data, including the same eight and 
one additional object. One of the objects in common be¬ 
tween these two searches, Kim 2, had in fact been previously 
found by |Kim et al. |2015| using data from the Stromlo 
MW Satellite Survey. In addition to that, Ki m fe Jerjen| 
(2015b) have discovered yet another object using Y1 Al data, 


Horologium II. The Panoramic Survey Telescope and Rapid 
Response System 1 (PAN-STARRS) has also been respon¬ 


sible for several recent discoveries of Milky-Way satellites 
(Laevens et al.|2014 Laevens fe Martin|2015 Laevens et al. 


20151. 


As described in |Becht ol et al. |2015|, several comple¬ 
mentary search strategies have been implemented within 
the DES Collaboration to search for stellar substructures. 
In |Bechtol et al.] ( |2015[ ), we used a conservative star selec¬ 
tion to ensure high stellar purity and completeness as well 
as a uniform held density over the survey footprint. The 
present work extends the results presented in |Bechtol et al.| 
(2015) by including stars at fainter magnitudes and by con¬ 


sidering a broader range of spatial extensions as well as ages 
and metallicites for the stellar populations composing new 
satellite systems. We also describe in detail the application 
of another search algorithm to the Y1A1 coadd data. To¬ 
gether, these analysis updates have enabled the discovery of 
a new candidate stellar cluster, DES J0034-4902, which we 
call DES 1, and the confirmation of Horologium II as a phys¬ 
ical stellar system. In Sj2]we describe the first-year DES data 
used. In Sj3]we describe the matched-filter algorithm applied 
to find the new systems. The new discovery is presented in 
)|4] In f|5]we report on the detection and characterization of 
Horologium II. Our final remarks are given in © 


2 DES DATA 

DES is a wide-held imaging program expected to cover 
about 5 000 deg 2 in the grizY bands down to ~ 24.6 th mag¬ 
nitude (at S/N ~ 10 for galaxies in g band; |The Dark En¬ 
ergy Survey Collaboration (2005)) in the southern equato¬ 


rial hemisphere for a period of hve years. It uses the Dark 
Energy Camera (DECam), a 3 deg 2 (2?2 diameter) mosaic 
camera with (//263 pixels on the prime focus of the CTIO 
Blanco 4 m telescope (Flaugher et al. 2015). The DECam im¬ 
ages are reduced by the DES Data Management (DESDM) 
team, which has developed a pipeline to process the data 
from basic single exposure instrumental corrections all the 
way to catalogue creation from calibrated coadded images. 
Here we use the DES year one coadd catalogue data (Y1A1), 
taken from August 2013 to February 2014. For more details 
on Y1A1 and DESDM we refer to |Sevilla et al.| ( |2011[ ), |Mohr| 
et al. (20121, and Gruendl et al., in prep. The stellar sample 


used in this work was drawn using the SExtractor param¬ 
eters FLAGS, SPREAD_M0DEL, and PSF magnitu des (|Bertin fe| 
Arnouts 1996 Desai et al. 2012 Bouy et al. 20131. These 


latter quantities are weighted average values taken from in¬ 
dividual exposures of each source. We used a source quality 
criterion of FLAGS ^ 3 over the gri hlters. 


As mentioned above, Kim & Jerjen (2015b I have discov 


ered one stellar object (Horologium II) in Y1A1 data that 
was not initially identified by Bechtol et al. (20151, or by 


Koposov et al. (2015 I. We believe that a primary reason for 


the non-detection of this object is that most of the probable 
member stars are fainter than g ~ 23 mag, which is where 


Bechtol et al. (20151 set the faint-end threshold to search 


of stellar objects. This conservative threshold to select stars 
was set to ensure high stellar purity and completeness, as 
well as a uniform held density over the survey footprint. 

In this work, we adopt a selection in SPREAD_M0DEL in¬ 
tended to increase stellar completeness, specifically, i-band 
jSPREAD_M0DEL| < 0.003 + SPREADERR_MODEL. A bright (faint) 
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g magnitude limit of MAG_PSF = 17 (MAG_PSF = 24) was also 
applied. The faint limit is 1 mag deeper than used by | Bechtel] 


et al. (2015). In order to prevent point sources with extreme 


colours (including red dwarfs from the Galactic disk) from 
contaminating the sample, a colour cut at —0.5 ^ g — r ^ 1.2 
was also used. 


Table 1. Parameter grid used to simulate SSPs for the search of 
star systems in DES Y1 footprint. 


Parameters 

Lower limit Upper limit 

steps 

log(Age) 

9.0 10.2 

0.3 

Distance (kpc) 

10 200 

10 

Metallicity, Z 

0.0002, 0.001 and 0.007 



3 SEARCH METHOD 


As discussed in Bechtol et al. (2015), several independent 
search methods were used in the original analysis of Y1A1 
data. In this section we describe in detail a different method, 
which was the one primarily used in this work. 


3.1 Matched Filter 


The Matched Filter (MF) technique has several applications 
for signal processing. In the context of astronomy, it has 
been used to detect low-density features and populations in 
imaging data (Szabo et al. 2011 Rockosi et al. 20021. We 
here use it to search for new star clusters and dwarf galaxies 
following on the work by Balbinot et al.| (2011). 

The number of stars as a function of position on the sky 
(a, 8) and of colour (c) and magnitude (m) may be generally 
described as 


n(a, 5,c, m ) = n c i{a, 5,c, m) + n bg (a,8,c,m). (1) 

The first term on the right hand side corresponds to the con¬ 
tribution by the cluster (cl) we want to discover, whereas the 
second term represents the background ( bg ), which includes 
foreground halo stars and background unresolved galaxies. 
We then split these terms into a normalization term and a 
probability distribution function (PDF): 

n c i(a,S,c,m) = ( C i(a,5)f c i(c,m ), (2) 

where £ c i and f c i are the number normalization and PDF 
on the color-magnitude diagram (CMD) plane, respectively, 
for the stellar population to be found. The stellar population 
may be extended in space (as in a stream), but we explicitly 
assume that its CMD is the same everywhere. As for the 
background stars, Galactic structure models show that both 
the number density and CMD vary as a function of position 
across the sky. So, we write 

n bg (a,5,c,m) = ( bg (a,6)f bg (a,5,c,m ). (3) 

With the definitions above, equation 0 then becomes 

n(a,5, c, m) = ( c i(a,5)f c i(c,m) + Cb g (a,5)f bg (a,5,c, m). 

(4) 

We bin stars into spatial pixels of area of 1(0 x l.'O, indexed 
by i, and color-magnitude bins of 0.01 mag x 0.05 mag, in¬ 
dexed by j. Details on the construction of the f c i and f bg 
PDFs are found in §3.2 and §3.3, respectively. With this 
notation, 

n(i,j) = Cd(i)fci(j) + C b g (i)f bg (i,j )■ (5) 

The left hand side is the expected number of stars in a given 
spatial pixel and CMD bin. If the actual number of stars 
observed in a catalog is N(i,j ), the variance between data 
and model is 

* W V MOMbi) • (6) 


The term in the denominator expresses the expected Pois¬ 
son fluctuation in the star counts, which, for simplicity, we 
assume to be dominated by the background. Minimizing the 
variance and solving for £ c ;(i), we have the number of ob¬ 
served stars that, according to the model given by equation 
Q, are consistent with the model. 

,_E j N (iJ)fci(j)/fb g (i,j) ( bg (i) 

CciW ~ Ej ) Ej f 2 Aj)/fb g (i,j ) ■ 

( 7 ) 

The output of the filter application is a stellar density map 
of stars which are probable cluster members stars, i.e., £ c i(i)- 
In practice, f bg (i,j) is generated from our target stellar cat¬ 
alog itself. We do that under the assumption that the con¬ 
tamination by any yet to be detected cluster, dwarf galaxy, 
or stellar stream, does not change the background PDF. As 
for the object PDF, we make use of simulated samples, as 
described in the next subsection. 


3.2 Model Grid 


Since we do not know a priori what stellar populations 
we will find, we create a grid of simple stellar popula¬ 
tions (SSPs) with the code GenCM!*] GenCMD uses PARSEC 


isochrones by Bressan et al. (2012) for different assumed 


distances and randomly selects stellar masses following a 
predefined initial mass function (IMF). Currently, we are 
adopting a Kroupa (2001) IMF for that purpose. Given each 


stellar mass, we interpolate among the isochrone entries to 
draw absolute magnitudes in the desired filters. These are 
converted into measured apparent magnitudes using the as¬ 
sumed model distance, magnitude uncertainties taken from 


Y1A1, and the reddening map of Schlegel et al. (1998). 


Positions on the sky may also be simulated assuming dif¬ 
ferent profile models. We simulate several SSPs at various 
ages, metallicities and distances covering a broader range 
of isochrone choices, including younger and higher metallic- 
ity stellar populations, than those adopted in |Bechtol et al.| 
( 2015| ). The parameter grid of these simulations is presented 
in Table Q] 


3.3 Object Detection 

We apply the MF method as presented in §3.1| to the stellar 
catalog using each of the SSPs in the model grid described 
in §3.2| In practice, the sky is partitioned into ~ 10° x 10° 
cells to account for local variations in the background CMD, 
which is empirically derived from the stars in each individual 
cell. This procedure generates one density map for every sky 
cell and for every point in the model grid. 

The maps are then convolved with Gaussian spatial ker¬ 
nels of different sizes {a = 0.'0 (no convolution) to a = 9.'0} 

1 https://github.com/balbinot/gencmd 
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Figure 1. DES coadd image cutout of DES 1 taken from the DES 
Science portal. The image is 1.'78 X 1.'78 centered on DES 1. The 
R,G,B channels correspond to the i, r, g bands. 


to highlight substructure on scales typical of star clusters 
and ultra-faint dwarf galaxies. In particular, smaller spatial 
kernels are suitable for the detection of more compact stel¬ 
lar systems. Our choice of spatial kernel sizes is such as to 
complement the ones adopted by Bechtol et al. (20151. 

As it is not practical to visually inspect all the resulting 
maps from the large number of combinations of sky cells, 
SSP models, and spatial convolution kernels, we use the 
SEXtractor code to automatically search for density peaks. 
In fact, the convolution kernels are applied from within 
SEXtractor itself as we run it on maps of different sky re¬ 
gions resulting from different SSP models. Any object found 
by SEXtractor in each map is recorded. Objects are then 
ranked according to the number of times they are detected 
by SEXtractor. This is done separately for each sky cell 
and convolution kernel. The SEXtractor parameters for the 
search were defined as those that maximized the recovery 
of simulated objects of different sizes and richness inserted 
into the DES stellar catalogue. 

The first ten objects in each region of the sky and for 
each convolution kernel are visually inspected to identify the 
most likely candidates. We visually checked the stellar den¬ 
sity map around them, the Poisson statistical significance 
above the background represented by their associated stel¬ 
lar peak, their number density profile and CMD. The density 
maps, significance and density profiles provide a basic assess¬ 
ment of the overdensities being found. The CMD helps us 
judge if this overdensity is consistent with a stellar popula¬ 
tion. All these diagnostic tools are shown in the next section 
for DES 1. 

We validated the detection method described above, 
which we call SparSEx, by applying it to simulated SSPs 
superposed on real SDSS and DES data. We also tried to 
recover faint MW satellites previously discovered in SDSS 
data. In particular, we have chosen seventeen objects found 
in SDSS data (eleven dwarf galaxies and six star clusters) 
that are characterized as distant and ultra faint objects. 


SparSEx detected all these stellar objects at the top of the 
object ranking lists. The rate of success for simulated stel¬ 
lar systems was the same. Finally, SparSEx detected all eight 
satellites reported by Bechtol et al. (20151 and Horologium II 
reported by Kim & Jerjen (2015b I. This later is discussed 
in more detail in |5| The ninth object detected by |Koposov| 


et al. (2015 I is in a region of Y1 data that is not included in 


the Y1A1 coadd that we searched. 


4 DES 1 


DES 1 stands out as the most conspicuous new candidate 
from our search. It is also directly seen as an overdensity of 
blue stellar sources in the DES coadd images (Figure [l]). In 
Figure [2] we show the number density of stars on the sky 
around this object (top panels). The left panel shows all 
classified stellar sources, as described in n and the mid¬ 
dle one shows only those close to the best-fit isochrone (see 
below). A clear overdensity is seen in both. The Poisson sig¬ 
nificance profile shows a very pronounced peak at about 1 (0 
from the object centre. It is built by taking the ratio of the 
number of stars internal to each radius r and in excess of 
the background (A7 &gd ), N a bj, to the expected fluctuation in 
the same background. N 0 bj = ( N 0 b s — Nb g d), where N 0 b s is 
the total number of observed stars. Nb g d is computed at an 
area ring at 30)0 < r < 34(0 from the center of DES 1. The 
bottom panels of Figure [2] show that there is not an over¬ 
density of galaxies at the position of DES 1, and therefore 
it is unlikely that misclassified galaxies can account for the 
apparent stellar overdensity. 

We use a maximum-likelihood technique to infer struc¬ 
tural and CMD parameters for DES 1. For simplicity, we 
assume that this object follows an exponential spatial pro¬ 
file with five free parameters: central RA (an) and DEC (do), 
position angle 9, ellipticity e, and exponential scale r e . We 
follow the same method as Martin et al. (|2008]) to find the 
best-fit solution and the parameter uncertainties. We then 
use the CMD of the most likely members of the system to fit 
an isochrone model, whose free parameters are: age, [Fe/H], 
(m — M) o, and Ay. The method is based on finding the peak 
likelihood in a series of model grids, as described in detail 
by Pieres et al (2015), in prep. 

In Figure [3j we show the results of the profile fit to 
DES 1. The first three panels show the likelihood values pro¬ 
jected on individual planes of this five-dimensional space, 
which all show well-defined peaks. The corresponding pa¬ 
rameter values and their uncertainties (computed as dis¬ 


cussed in Martin et al. (2008) and Pieres et al 2015, in 
prep.) are listed in Table^2pThe last panel shows the in¬ 
dividual stars coded by their membership probabilities. We 
note that DES 1 is a quite elongated object (e ~ 0.7). Fig¬ 
ure [4] shows a binned radial density profile compared to the 
best fit exponential model. The central density of DES 1 is 
~ 200 stars/arcmin' 2 . A clear excess of stars relative to the 
background is seen out to ~ 2(0. In fact, a systematic ex¬ 
cess of stars relative the best-fit model is seen beyond ~ 1(0, 
which may indicate that a King profile may better describe 
the observed profile. 

The CMD for this candidate is shown on the left panel 
of Figure [5] Only stars inside the radius corresponding to 
the peak in the Poisson significance profile shown in Figure 
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Figure 2. Top left panel: on-sky number density map of stellar sources around candidate DES 1. All stars are included. Top middle panel: 
similar to previous panel, but now only stars which lie close to the best fitting isochrone shown in Figure [5] are used. Top right panel: 
Significance as a function of radial distance from the centre of DES 1. The solid (dashed) line correspond to all (isochrone filtered) stars 
as indicated. The corresponding panels at the bottom show the same plots but using the distribution of sources classified as galaxies. 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

log(r/arcmin) 


Figure 4. Solid points show a binned version of the radial density 
profile of DES 1. The error bars are 1-sigma Poisson. The dashed 
line is the best maximum-likelihood fit model to an exponential 
profile. 


[2]are shown. The middle panel shows the field CMD in a ring 
of equal area, whose inner radius equals ten times the radius 
of the object’s CMD. The best-fit isochrone is also shown; 


parameter values are listed in Table [2] Also shown are the 
sequences bracketing the best-fit isochrone fit at a distance 
of VO.l 2 + MAG_ERR 2 + C0L_ERR 2 on the CMD plane, where 
MAG_ERR and C0L_ERR are the mean photometric uncertain¬ 
ties along the CMD axes. We use the CMD space between 
them to filter the most likely cluster stars. See the middle 
and right panels of Figure [2] for a density map and a Poisson 
significance profile of objects inside this isochrone filter, re¬ 
spectively,. The CMD difference relative to the background 
field is shown as the Hess diagram in the right panel of Fig¬ 
ure [5] The main sequence turn-off (MSTO) and sub-giant 
branch (SGB) are clearly visible. 

As previously mentioned, we summarize the inferred 
properties of DES 1 in Table [2] The table lists positions, 
structural parameters, central (cr c ) and background (cr bs ) 
densities, half-light radius ( rh ), distance {Dq), Absolute 
Magnitude {My), Test Statistic (TS), and peak Poisson sig¬ 
nificance (PS) value, as well as best-fit CMD parameters. 
The TS is based on the likelihood ratio between a hypoth¬ 
esis that includes an object versus a field-only hypothesis 
(see [Bechtol et al.|2015| ). 

A maximum-likelihood approach that simultaneously 
fits the profile (assuming a Plummer model) and the dis¬ 
tance (but assuming an age of 12 Gyrs and a spread in metal- 
licities) has also been tried. It yields a distance modulus of 
(m — M )0 = 19.6, in agreement with the method described 
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Figure 3. Upper left panel: likelihood map for DES 1 projected onto the position angle and eccentricity plane. Upper right: likelihood 
map for DES 1 projected onto the position angle and exponential scale plane. Lower left: likelihood map for DES 1 projected onto the 
central equatorial coordinates plane. Lower right: Spatial map of stars with probability larger than 1% to belong to DES 1 color-coded 
by probability. The best-fit parameters are listed in Table [2] 
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Background 
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Figure 5. Left panel: CMD of stars within 2/0 from the centre of DES1. The maximum-likelihood isochrone fit is shown, along with 
ridge lines meant to bracket the most likely members. Middle panel: CMD of background stars in an equal area on the sky as the previous 
panel. Right panel: Hess diagram of the CMD difference between the two previous panels. 


earlier. The alternative rh, however, is larger, rn — 1?0. 
Visual fits were also independently made to the object’s 
CMD. Again, the results agree well with those from the 
maximum-likelihood fits shown on the table: log(Age) = 9.9, 
[Fe/H] = -1.98, A v = 0.03, and [m - M) 0 = 19.9. The 
quoted My estimate was computed by integrating over all 
masses along the best fit model isochrone assuming a |Kroupaj 
(2001) IMF, and normalizing the number of objects by those 
observed in the CMD with r < 23 mag and which fall in the 
isochrone filter. We convert from DES g and r magnitudes 


to V magnitudes using 

Ssdss = Sdes + 0.104( 5d es — ?~des) — 0.01 
fSDSS = tdes + 0.102(<?des — »~des) — 0.02 (8) 

V = psdss — 0.59(psdss — Usdss) — 0.01. 

These transformations were derived by |Jester et al. ([2005]) 
using a SDSS stellar calibration sample. 
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Figure 6. All panels are the same as those in Fig. [2] but now for the Florologium II satellite dwarf. 


Table 2. Properties of DES 1 


Parameters 

DES 1 

ao (J2000) 

<5 0 (J2000) 

r e 

9 (deg) 
e 

ct c (stars/arcmin 2 ) 

(T bg (stars/arcmin 2 ) 
rh 

Dq (kpc) 
rh (pc) 

My 

TS 

PS 

00 h 33 m 59.7 s ± 1.7 s 
—49°02'20.0" ± 13.7" 

Of 23 ± O.'IO 
-57.9 ±22.2 

0.69 ±0.22 

204.98 ± 35.81 

2.02 ±0.01 

Of39 ± 0.'17 
~ 87.1 

9.88 ±4.31 

_Q qc+0.69 
o.UO_ 0 42 

134.7 

11.7 ±3.1 

[Fe/H] 

-1.98 

l°g (Age) 

10.00 

Ay 

0.0 

O 

S' 

J. 

19.70 


5 HOROLOGIUM II 


As mentioned earlier, Kim & Jerjen (2015b) report on the 
identification of an extra MW satellite galaxy candidate be¬ 
sides those found by |Bechtol et al.| ( |2015| ) and |Koposov et ah] 
(20151. Our reanalysis of Y1A1 presented here confirms this 
object, Horologium II, as a real stellar system. In fact, once 
we allow for a deeper magnitude threshold, as explained in 
fj2] we detect it not only with the method described in Sj3]but 


also with the maximum-likelihood satellite search method 
detailed out in jBechtol et al.| ( |2015| ), 

Figure [6] shows the same information as Figure [2] but 
now for Horologium II. A clear overdensity of stars is seen 
in the density map on the sky. The statistical significance 
of this overdensity is close to 8 times the expected Poisson 
fluctuation in the background. 

Using the same maximum-likelihood fit method to the 
on-sky and CMD distributions as in the previous section, we 
obtain the following structural and isochrone parameters: 
(m — M )o ~ 19.5 (heliocentric distance D ~ 79.4 kpc), 
log(Age) ~ 9.9, [Fe/H] ~ -1.2, A v ~ 0.01, r h ~ 2.2' (51.2 
pc), position angle 0 ~ 106°, elipticity e — 0.6, and absolute 
magnitude My ~ —3.1. 


Our distance, size and flatenning estimates are in very 
good agreement (within la) with those from the discovery 
paper by Kim & Jerjen ([2015b). The position angle and 
absolute luminosity are within 2a of their quoted values. 
The largest discrepancies occur for the isocrhone parame¬ 
ters. We fit the observed CMD of Horologium II to a younger 
and more metal rich PARSEC model than Ki m fc Jerjenj 
( |2015b| ). However, our likelihood function over the metallic- 
ity and age plane exhibits a tail of high likelihood values 
towards lower Z and older ages. The discrepancy may also 
be related to the fact that our listed redenning value comes 
out of the maximum-likelihood CMD fit, whereas the val¬ 
ues from Schlegel et al. (1998) maps {with corrections from 
|Schlafly~fc Finkbeiner] (2011)} are used in the discovery pa¬ 
per. 
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6 CONCLUSIONS 


In this paper we make a deeper probe on the DES Y1A1 
catalog in search for additional Galactic satellites besides 


those previously reported by the collaboration (Bechtol et al. 


20151. We report the discovery of a new stellar system in the 


MW halo, using catalogs based on first-year data from the 
Dark Energy Survey. We have explored the data at least 


1 mag deeper (g < 24 mag) than previously done in Bech 


tol et al. (2015). The candidate adds to the ten previously 


identified systems found using DECam images by |Bechtol 


et al. (2015); Koposov et al. (2015); Kim et al. (20151; 


Kim 


& Jerjen (2015b I. We also confirm the dwarf galaxy candi¬ 


date Horologium II, originally discovered by |Kim fc Jerjen| 
12015b), as a significant overdensity in the Y1A1 catalog. 


Our best fit structural parameters to this later are in gen¬ 
eral agreement with the ones derived by those authors, al¬ 
though the isochrone fit tends to point to a younger and 
more metal-rich object than previously reported. 

DES 1 is detected as a significant stellar overdensity 
both spatially and on the CMD plane. Isochrone fits based 
on two very different methods show that it is made up of 
old and metal poor stars, as commonly observed in MW 
satellites found in the Galactic halo. Extinction towards the 


object, as estimated from the Schlegel et al. (1998) maps, is 


found to be relatively small (~ 0.03 mag). 

Maximum-likelihood fits of the spatial profile of DES 1 
yield a scale radius of r e ~ 0(2, which at a distance ~ 87kpc 
corresponds to a physical size of rh — 9.88 pc. Its estimated 
distance places this faint cluster candidate as one of farthest 
away from the Sun. The absolute magnitude has been de¬ 


termined using a similar approach as Koposov et al. (20151, 


yielding My ~ —3.05 for DES 1. Taken together, the phys¬ 
ical size and luminosity place DES 1 in the locus occupied 
by low-luminosity star clusters. DES 1 is also significantly 
elongated (e ~ 0.7), something that is apparent not only 
from the profile fit, but also from the distribution of the 
stars on the sky (Figures [I] and [3]). It is, in fact, the most 
elongated halo cluster known to date, although, given the er- 
rorbars, its excentricity is marginally consistent with those 
of Kim 1, Laevens 3 (see Kim & Jerjen|2015a Laevens et al. 


2015| respectively). The very high inferred excentricity sug 


gests that DES 1 is in dynamical process of tidal disruption, 
despite its large distance, and makes it a very interesting 
object for deeper imaging and spectroscopic follow up. 

Due to the low number of probable member stars de¬ 
tected in the DES imaging, it is difficult to extract more 
reliable information about DES 1 at this stage. Steps to ac¬ 
quire deeper imaging of this object are already under way. 

A search for satellites in data collected by the Dark 
Energy Survey during its second season, including new areas 
sky, is under way. It is likely that additional new stellar 
systems will be discovered soon. 
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